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INTRODUCTION 



The original planning and experimentation from which these programs 
evolved was done in a Math/Physics seminar at Hanover High School in the 
Spring of 1969. We wished to explore and understand the physical forces 
involved in orbital flight at a depth for which no secondary school literature 
appeared to exist. The DTSS* presented an interesting opportunity for the 
students to do actual constructive research using mathematical models and 
computer simulations. 

The mathematics is essentially that of a superior Algebra II course that 
stresses Trigonometry and Analytic Geometry, A student's ability in these 
areas may be sorely tried and perhaps extended if these programs and supporting 
material are dissected and examined. Many problem solutions normally ob- 
tainable only through the methods of calculus are found through computer 
programs utilizing some interesting algorithms. The original questions asked 
by the seminar students were non-trivial, and the resulting analyses were, 
consequently, slightly more sophisticated than was expected. Despite this 
creeping complexity, the questions that were answered at all were answered 
with a substantial degree of integrity. 

The physics represents an extension of the basic physics employed in 
the excellent Holt, Rinehart & Winston paperback series on space flight with 
specific emphasis on "Mathematics of Space Flight", the most utilitarian of 
these paperbacks. Obviously, the programs go beyond this introductory 
material but may still be considered a linear extension of it. Vectors and 
energy levels are the prime tools of kinetics analysis • At least one seminar 
of a good physics course such as PSSC is heartily recommended as a pre- ' 
requisite to this material although much of the necessary background is well 
introduced by the paperback text. 

The forces and energies studied have such a large magnitude and many 
of the changes sought after are so minute that the computer cannot give the 
accuracy of computation that some students would desire. However, the 
numerical analyses yield data that is sufficient for most uses. 

During the Spring of 1970, some Hanover students are planning to write 
exercise sets around the present programs, and perhaps, some new programs 
as well. These exercises will lie in three categories: (1) questions requiring 
extensive knowledge of the given program and its accompanying Topic Outline 
and also requiring many directed runs of the present program; (2) questions 



* Dartmouth Time-Sharing System 



requiring minor re-writing of the program to present data from a different 
point of view; and (3) compilation of tables to aid in analysis directed toward 
writing new exploratory programs and requiring continuous and systematic 
re-programming. These exercise sets may be available in a later edition of 
this topic outline. The following program notes may help in creating ad hoc 
exercises. 

(1) Elliptic Orbits : This topic is a time analysis of closed orbits. 
The time analysis is made possible through the use of a computer 
algorithm that rather accurately replaces the commonly used 
elliptic integrals. The laws of planetary motion may be studied 
quite closely after minor substitutions are made in the program. 
Keppler's laws, as one example, are directly observable on a 
time plot. 

(2) Rendezvous : The simplest form of the two ship rendezvous is solved 
and the time function is used again to track the interception in 
either tabular or graphic displays. The program is stripped to the 
absolute minimum to emphasize "window" and relative orbital 
constants. Obviously, a reversal of the ship roles in the program 
leads to the Hohman transfer. The program as it stands gives a 
clear picture of circular and elliptic orbits of the same magnitude. 

(3) Orbital Transfer: This program computers the effect of adding any 
energy vector to any orbit configuration at any point. The potential 
of this prograni as yet is largely unexplored and its limitations and 
weaknesses have not yet been fully identified. Suggested exercises 
would include studies of orbital stability, relative influences of 

the multiple variables, and of course, problems in rotation of 
orbital axes. The advanced student may use an easy variation of 
the program as a flight simulator in exchanging orbits. 



Future developments will hopefully extend Rendezvous to include all manner 
of orbits. This will, of course, require the most careful exploitation of 
Orbital Transfer to control rotation and alignment of orbital axes. 
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I. ELUPTIC ORBITS 



REF: ELORB 



Time analysis of closed orbits has been a subject of fascination to astro- 
nomers, physicists, and mathematicians alike for some 300 years. Theoretically, 
if the orbit is known, then the orbiting body's position is strictly a function of 
time. Unfortunately for the untutored amateur or elementary student of orbital 
mechanics this position plot as a function of time gives rise to some complicated 
applied mathematics. 

The advent of the high-speed computer offers an alternative to the elliptic 
calculus developed by the giants of the 18th and 19th centuries to handle 
"celestial mechanics. The ability of a computer to handle large numbers of 
complex calculations in fractions of a second allows "back solving" of systems 
of equations which, if solved in a straightforward manner, lie in the domain of 
advanced mathematics. 

Before examination of the actual equations from which the program is con- 
structed, a look at the physical situation is convenient. In Figure 1 is shown 
a typical elliptic orbit with Earth at its principle focus. A space ship is shown 
at an instant of time along with the area which its radius vector has swept out. 
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Since an orbit which is undisturbed by additions or subtractions of energy 
has constants for its descriptive parameters, it is obvious that a few astronomical 
observations will yield: 

1. Closest approach to Earth 

2. Furthest recession from Earth 

From these data all other significant parameters may be calculated. From an 
arbitrary zero time it seems as though after K seconds a given area is swept 
out by the radius vector which is directly dependent on 9, the radius vector 
angle. 

Keppler established that the radius vector sweeps out equal areas in equal 
times. If it were easily calculable, the area of an ellipse in terms of 9 would 
be of immense aid since we could then directly relate t and 9. In Figure 2 
such a calculation is established. The result is an integral which is included 
in the text as a "given" equation but should be briefly examined at this point. 




Figure 2. The Differential Triangular Area of an Ellipse 



From Figure 2 the area of A AOB is ir.r.d9, and the integral (where r = L/l+e cos9) 
is 

d9 

a (1 + e cos9)2 



2 / < 



which is a standard transcendental with a solution readily available from any 
set of integral tables. 
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In most mathematical texts the solution to the area problem is mani- 
pulated so that t = f(e) which is, of course, the alDsolute inverse to what 
is required, i.e. , 9 = f(t). It is precisely at this point that a secondary 
school student with a computer can match the calculations of the more ad- 
vanced student. 



The program associated with the mathematical portion of this topic 
outline can be set to report angular position (r,9) given any time interval. This 
is accomplished by having the computer search for a 9 that will produce the 
required time. The present. program is adjusted for 5 minute intervals with 
an accuracy of .001 minutes, or 15 search sweeps, whichever occurs ilrst. 
Tnis avoids the inverse solution to the differential triangular area of the 
ellipse which leads to elliptic integrals which would haye to be solved for 
each time plot. . 

When "approach" and "recession" distances are given to the computer in 
response to its input interrogation, the program not only presents most of the 
orbital data in tabular form but on command will track the orbiting object giving 
position plots with velocities at 5 minute time intervals'. 

As a consequence, any student with an elementary background in mathe- 
matics who gain access to a time-shared computer, may study closed orbits 
with their associated laws of physical behavior. The more serious student who 
cares to pursue the mathematics of this physical behavior will find all the 
equations necessary for further exploration and possible re-programming. 

Objects in free fall around Earth fall in orbits which are mathematically 
describably as conic sections. The principal focus is at Earth center. The 
polar coordinate system is used to describe a ship's orbital position. The 
reference axis is through Earth center on the orbital perigee. The ship's 
position in an elliptic (as well as any other conic orbit) is then 



1 + e cos 9 



(El) 
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Figure 3, Orbital Parameters (Elliptic Orbit) 



A: Apogee 
P: Perigee 

L: Semi Latus Rectum 



r: Radius vector 
9: Radius angle 
e: Eccentricity 



a: Semi Major Axis 
b: Semi Minor Axis 
ae: Semi Focal Length 

NOTE: Counterclockwise orbit is with Earth rotation and will be standard. 
The distance parameters in terms ot A, P are: 



a 



(A + P)/2 



(E2) 



e 



(A - P)/(A + P) 



(E3) 



L 



a(l - e2) 



(E4) 



b 



a >i 1 - e2 
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The inertial period, Pe, of an elliptic orbit is 
Pe = 



V Gm 



seconds 



where Gm, the Earth gravitational constant is 

Gm = 62747 nautical miles/second^ . 

The ship's velocity at any point in orbit is 
V = ^ Gm{2 



2.) nmps, 
a 



Velocities at apogee and perigee will be 
vp 



. / Gm ( 1 + e ) nmps 
V A 



v^ = sjom (1 - e ) 



nmps . 

And the area swept out by the radius vector is 

e sin 9 + 2 



A(r,9) = L?. 

2 



(e^-Dd+e cos 9) (1 - e2)3/2 

By Keppler's 2nd Law, the following relationship holds 

A(r, 9) = t seconds 
k 

If (Ell) is solved for one complete revolution 

A(r,9) = /Tab 
and the following equalities may be established: 

Tfab = 2Tr a^ k 



Tan"l ( ^l-e2 Tan 9) 

2 



Gm 



aGm 



60bGm 



e cos 9 + 2 Tan"-^ (v/l-e^ Tan9) 

(e^-l)(l+e cos 9) (l-e2)3/2 2 



(E6) 



(E7) 



(E8) 



(E9) 
(ElO) 

(Ell) 



(El 2) 



(E13) 



(E14) 



minutes 

(El 5) 



To directly solve (E15) for 9 in terms of t_ is evidently quite difficult, so the 
time function is computed through a computer convergence or search program. 
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ELi3HB 03/ 18/70 09: 51 



THIS PRCIGKAM IS KEYED T0 THE T0PIC 0UTLINE 'ELLIPTIC 
aKBITS' 1HN0UGH THF. REM C0HJMN REFERENCES 

CLOSEST APPR3ACH T0 EARTH? 150 
FURTHEST RECESS! 3N FROM EARTH? 9 50 



03NSTANT PARAMETERS 



SEMI MAJOR AXIS IS 3992 
SEMI MINOR AXIS IS 397 1.9 1 
SEMI LATUS RECTUM IS 3951 .92 
SEMI FOCAL LENGTH IS 400. 
ECCENTRICITY IS 0.1002 
INERTIAL PERIOD IS 105.443 
PERIGEE VELOCITY IS 1578 2.2 
APOGEE VELOCITY IS 12907.5 



VARIABLE PARAMETERS 



ELAPSED TIME RADIUS ANGLE RADIUS V VELOCITY 



0 


1.09863 E-2 


3592. 


1 5782. 2 


5 


18 .3 525 


3609 . 64 


15711.9 


10 


37. 5403 


3661 . 04 


1 5 509. 3 


1 5 


55.9314 


3741 .88 


15196. 7 


20 


73.905 


38 45. 1 1 


1 4807.8 


25 


9 1 . 3733 


3961 . 43 


1 4382. 3 


30 


108. 29 2 


4080. 24 


1 39 60. 5 


35 


124. 706 


419 1 . 02 


1 3573. 


40 


1 40. 68 


428 4. 


1 3264. 2 


45 


1 56. 28 1 


4351 . 07 


1 3041 . 7 


50 


17 1.661 


4386.84 


1 2924. 3 


55 


186.976 


4388. 38 


129 19 . 3 


60 


202. 346 


4355. 58 


1 3026.9 


65 


21 7.9 36 


429 1 . 03 


1 3240. 7 


70 


233.855 


4200. 1 5 


1 3546.8 


75 


250. 225 


4090. 6 


1 39 2 4. 4 


80 


267. 1 1 1 


3971.98 


1 4344. 4 


85 


28 4. 535 


3854.98 


1 4771 . 2 


90 


302. 465 


3750. 21 


1 51 64.9 


95 


320.8 12 


3667. 1 2 


1 548 5. 5 


100 


339 • 489 


3612.86 


1 5699 . 1 


105 


356 • 33 


3592. 1 4 


1 578 1 - 6 



TIME: 2.000 SEC* 
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100 PRINT "THIS PROGRAM IS KEYED T0 THE TOPIC OUTLINE 'ELLIPTIC" 
110 PRINT "ORBITS' THROUGH THE REM COLUMN REFERENCES" 
120 PRINT 

130 PRINT "CLOSEST APPROACH TO EARTH"; 
lAO INPUT X 

150 LET X = X+3442 'P 
160 PRINT "FURTHEST RECESSION FROM EARTH"; 
170 INPUT Y 

180 LET Y = Y+3442 'A 
190 LET A = (X + Y)/2 •CE2) 
200 LET E=(Y-X)/(Y+X) 
202 IF E< 1 THEN 210 

204 PRINT "ECCENTRICITY OF SELECTED ORBIT I S**E 
206 PRINT "REDUCE ECCENTRICITY BELOW E= 1 " 
208 GOTO 120 

210 LET L = A*(l-Et2) 'CE^) 
220 LET B = A*SQR(L/A) •( E5) 

230 LET P = 2*3« 1 41 59265*SQRC At 3/62747) •CE6) 
240 LET VC2)= SQRC 627 47/Y ) * SQRC 1 - E) 
250 LET VC 1) = SQRC 62747/X)*SQRC 1+E5 
260 PRINT 
270 PRINT 

2B0 PRINT "CONSTANT PARAMETERS" 

290 PRINT "- " 

300 PRINT 

310 PRINT "SEMI MAJOR AXIS I S"A 

320 PRINT "SEMI MINOR AXIS I S"B 

330 PRINT "SEMI LATUS RECTUM I S"L 

340 PRINT" SEMI FOCAL LENGTH I S"A*E 

350 PRINT"ECCENTRI CI TY I S"E 

360 PRINT "INERTIAL PERIOD IS" P/60 

370 PRINT "PERIGEE VELOCITY IS" V(l)*3600 

380 PRINT "APOGEE VELOCITY IS" V(2)*3600 

390 PRINT 

400 PRINT 

410 PRINT "VARIABLE PARAMETERS" 

420 PRINT " " 

430 PRINT 

432 PRINT "ELAPSED TIME"> "RADI US AN GL E"> "RADI US V"> "VELO CI TY " 

433 PRINT " ">•• 

435 PRINT 



440 


LET 


N ' 


INT(P/60) 




450 


FOR 


M = 


0 T0 N STEP 5 


'SEARCH L00P 


455 


LET 


3=0 






460 


LET 


0 = 


2*3. 1 41 59 26 5 


•0VER 


470 


LET 


U = 


0 


•UNDER 


4B0 


LET 


z = 


(0+U)/2 


•TRIAL ANGLE 


490 


LET 


T( 1) 


= Lt 2*SQK( A*62747)/( 60*B*62747) 


•( El 5) 


500 


LET 


T( 2) 


= E*SINCZ)/C (Et2- 1)*( 1+E*C0SCZ) ) ) 


•( El 5) 


510 


LET 


T( 3) 


= 2/SQR( 1-Et2>t 3 


•( El 5) 








is 
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( C0NTIN5UED) 



520 LET T(4) = ATN( SQR( 1-Et 2)*TAN(Z/2) ) '(£15) 
525 LET S=S+ 1 

530 LET T = T( 1)*(T(2) + (T(3)*T( 4))) '(ElS) 

532 IF Z <= 3. 14159265 THEN 540 

535 LET T=T+2*3. 1 41 59265*SQR( At 3/62747)/60 

540 IF ABS(T-M) < .001 THEN 600 'EXIT 
545 IF S=15 THEN 600 
550 IF T > M THEN 580 

560 LET U = Z 'NEWUNDER 
570 G0T0 48 0 

580 LET 0 = Z 'NEW 0VER 

590 G0T0 480 

600 LET R=L/( 1 + E*C0S(Z) ) 

610 LET V = SQR ( ( 2*A* 62747-R*62747) / ( A* R) ) ' VEL . 

„620 PRINT M*Z* 180/3. 1 41 59265* R* V*3600 

630 NEXT M 'END L00P 

700 END 
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IL RENDEZVOUS 

REF: CHASE 



The rendezvous problem has been simplified to the minimum necessary to 
illustrate problems of pursuit. Essentially different orbits have different velo- 
cities, and movement between them requires careful calculation if interception 
is the goal. The following figures illustrate the lead problem and its theoretical 
solution. 

Chase Ship (1) 

Descent 
Orbit 




Figure 1 . 



If the descent is from a co-radial position, the chase ship will be behind 
the target due to the difference in velocities and times betv^een the descent 
orbit and the target orbit. 
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Chase Ship (1) 



Lead Interval 
(Window) 



Figure 2 . 



The chase ship initiates its maneuvers froi.. a calculated lead interval 
based on co-solving the three orbits for elapsed time and average velocities. 



The descent orbit is interesting in its velocity changes. Initially, the 
chase ship sheds velocity as the elliptic orbit, with its apogee at the chase 
ship orbit, has a lower energy level at that point. As the ship descends with 
the force of gravity to its perigee at the target orbit, its velocity increases 
to the necessary energy level to "crack the whip" at perigee. So once again, 
velocity must be reduced. In essence, velocity is twice reduced, and in 
nearly the same amounts, to enter a faster orbit from a slower orbit. 

When the program is initiated the original ship orbits are randomized within 
certain bounds before being presented to the student. A Firing Table is generated 
which presents all the necessary'' navigational and maneuvering data for a ren- 
dezvous. If the student wishes to track the rendezvous, data is presented in 
intervals of the student's choice. The search program from ELORB is incorporated 
into the descent orbit computations and certain other data is computed from this 
to show the following: 
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1. Elapsed time in minutes 



2. Relative distance between ships in nautical miles 

3. Relative angle between the ships (in degrees) 

4. Relative velocity between the ships in nautical miles per hour 

5. Time in minutes until rendezvous 

One of the objectives of this program is for the student to become familiar 
enough with the program - physics and mathematics - to reverse the positions 
of the two ships. 

Suppose two ships are in differing circular orbits and a rendezvous is 
desired. Suppose the target ship is in orbit of radius r2 and the chase ship 
is in orbit of radius ri and further more, the chase ship is displaced L radians 
from the target and is in the larger (and slower) orbit. The chase ship must 
now establish a transfer orbit with its apogee on the chase ship original orbit 
and the perigee on the target ship's orbit. 

To describe the transfer orbit we have 



A = ri (El) 




Figure 3 . 



The velocity of the transfer orbit at apogee is 



/cm (1 

V n 



Vq = /Gm (1 - e) (E4) 



Hence, the amount of velocity to be lost by retrofire is the difference 
between the velocity of the chase ship's orbit and (E4) . 



XI = / cm - / Gm( 1 - e) 
\/ri V ri 



(E5) 
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The time to perigee in the transfer orbit is half the period of the new orbit 
or 



^ Gm 



At perigee the chase ship has a velocity higher than its tangential circular 
orbit and must again reduce velocity to match the target orbit. 



^2 



Gm(l + e) - / Gm 



If, for the purposes of simplicity, only the transfer time is considered, the 
chase ship must fire its motors precisely t-^ seconds before the target ship 
reaches the anticipated perigee rendezvous. 



Since the target ship travels 2 radians in 

seconds , 



t = 27r, / ^2^ 



Gm 



then its angular velocity may be considered 



rad/sec. 



r2^ 



Gm 



(E6) 



(E7) 



(E8) 



(E9) 



Hence, it must travel w^. radians to rendezvous 



w. 



^ = (E6)(E9) = TT./si 



(ElO) 



The chase ship must be radians from the perigee rendezvous so the 
chase ship's lead at time of firing must be 



w = Ti 



- 1 



^2' 



radians 



The relative angular velocity of the ships is 
V = 1 - 1 



Gm 



f r i3 
Gm 



rad/sec , 



(E12) 
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The amount of delay necessary for the relative ship velocities to change the 
relative ship angle to the proper ''window" is found from the two possible 
cases in Figure 4. 



Chase ship 

Window 




Tcirget 
ship 



Chase ship 

Window 




Case I 



Case II 



Figure 4 . 

In Case I the target ship is outside the window and 

t2 = (2 - W2 - w)/v . 
In Case II the target ship is inside the window and 

t2 = \{2 IT -W2) + (2^ -w)/v. 



The argument for a rendezvous with the chase ship occupying the inner 
position is similar. 



(E13) 



(E14) 
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(HASE 03/18/70 09:35 



FIRING TABLE 



CHASE SHIP ANGLE ISO 
(HASE SHIP RADIUS IS 49 9 3i67 
TARGET SHIP ANGLE IS 254.557 
TARGET SHIP RADIUS IS 4088.47 

WINDOW ANGLE IS 30.7019 
TIME T0 WIND0W IS 87. 5434 
REDUCE VELOCITY BY 652.627 
DESCENT TIME IS 63.9645 
REDUCE VELOCITY BY 686. 131 



D0 YOU WISH T0 TRACK THE SKIPS? YES 
TRACKING INTERVAL? 10 



TRACKING DATA 



ELAPSED 


RELATI VE 


RELATI VE 


TIME 


AN CLE 


DI STANCE 


0 


105. 443 


7247. 42 


10 


96.9051 


6823. 6 


20 


88. 3676 


6363. 1 


30 


79.8 3 


58 68. 77 


40 


71.2924 


5343. 75 


50 


62. 7549 


479 1 . 56 


60 


54. 2173 


4216.25 


70 


45. 6797 


3622. 66 


80 


37. 1 422 


3017*08 


87. 5A3A 


30. 7019 


2557.88 


EXECUTE INITIAL RETR0FIRE 




92. 5A3A 


26.888 


2279.87 


97. 5434 


23. 1 1 79 


1992. 76 


102. 543 


19 . 5237 


1707. 75 


107. 543 


1 6 . 1 49 2 


1 430. 23 


1 12. 543 


1 3. 0823 


1 1 68 . ^^7 


1 17. 543 


10. 367 


927.94 


122. 543 


7.95924 


709. 123 


127. 543 


5.9 4703 


520. 31 5 


1 32. 543 


4. 28 638 


362.216 


1 37. 543 


2.88941 


232. 557 


1 42. 543 


1 . 7561 


1 33. 1 47 


1 47. 543 


0. 75464 


54. 5755 


EXECUTE FliNiAL 


KETR0FI RE 




151. 50b 


2. 19 522 E-2 


1 . 56645 



RELATIVE 
VEL0CI TY 

1 342. 1 
1 342. 1 
1342. 1 
1 342. 1 
1 342. 1 
1342. 1 
1 342. 1 
1 342. 1 
1 342. 1 
1 342. 1 



19 58 . 61 
18 52. 1 A 
1 6 79 . 49 
J 448 • 53 
1 1 69 . 55 
856. 572 
528. 478 
203. 01 4 
98. 1257 
353. 32 
545. 032 
6 58 . 10 3 



2. 1 4577 E- 4 



TIME T0 
RENDE2 V0US 

151. 508 
141. 508 
131. 508 
121. 508 
111. 508 
101. 508 
9 1 . 508 
8 1 . 508 
71 . 508 
63.9645 



58.9645 
53*9 645 
-^.9645 
43^9 645 
38.9645 
33. 9645 
28.9645 
23.9645 
58.9645 
13.9645 
8.96452 
3.96452 



5.08 502 E-7 



TIME: 1.651 SEC 
HEADY 
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^0 



CHASE 



100 
105 
110 

120 
130 
140 
150 
160 
17C 
180 
19 0 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 

36 5 

37 0 
380 
39 0 
400 
410 
420 
430 
440 
450 
460 
470 
«0 
-^0 
500 
510 
520 
530 
540 
550 
560 
570 



RANDOMIZE 

LET J = 180/3. 141 59265 
LET R<1) = RND* 1000+4000 
LET R<-.) = RND* 1000+4000 
LET W< 1) = RND*6. 2831853 
LET U<2) = HND*6. 2831853 
LET V<1> = SQR< 62747/R( 1) ) 
LET V<2) = SQH< 62747/R<2) ) 
IF H< i) > RC 2) THEN 200 
IF R< 1) < RC2) THEN 100 
G0T0 100 
LET A = K< 1) 
LET P = RC 2) 
LET E = <A-P)/<A+P) 
LET S = " < A+ P) / 2 

LET TCI) = 3. 1 41 59265* SQR< St 3/62747) 

LET Vi = 3. 14159265*SQR(St3/R<2)t3)-3. 14159265 

LET V<3)=1/SQR(R( l)t 3/ 62747) 

LET V< 4)= 1/SQK( R<2) t 3/62747) 

LET V = V<4)-V<3) 

LET Z = 6.2831853-W< 1) 

LET W< 1) = 0 

LET U<2) = U<2)+Z 

IF V<C2) < 6.2831853 THEN 340 

LET W<2) = W(2)-6. 2831853 

IF U<2) < <6.28318530-V;) THEN 380 

LET Y = <6.28318530-VJ(2)) + <6.28318530-W) 

LET T<2) = Y/V 

IF T<2)/60>100 THEN 100 

G0T0 400 

LET Y = ( 6. 28318 530- W< 2) ) - W 
G0T0 360 

LET X<1) = SQR<62747/A)-SQR< 62747*< l' E)/A) 
LET X(2) = SQR<62747*< 1 + E)/P)-SQR< 62747/P) 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 



•0RI GINATE 



GATE 



•APOGEE 
•PERI GEE 
•ECCENT 
•SEMI MAJ 
•COAST TIME 
•WINDOW 



•REL 

•INI TI ALIZE 
•0RBI TS 



•TARGET ADV. 
•TIME TO FIRE 

•EXIT 

•TARGET ADV. 



•EXCESS 



'VEL0CI TIES 



•FIRING TABLE^' 



••CHASE SHIP ANGLE I S^^UC 1 ) * 18 0/ 3. 1 4 1 59 265 
"CHASE SHIP RADIUS IS"R<1) 

••TARGET SHIP ANGLE I S^' W< 2) * 18 0/ 3. 1 4 1 59 26 5 
••TARGET SHIP RADIUS I S'^RC 2) 

••WINDOW ANGLE I S^^ W* 13 0/ 3. 1 41 59 265 
••TIME T0 WINDOW IS'^T<2)/60 
••REDUCE VELOCITY BY"X< 1)*3600 
"DESCENT TIME IS^'T<l)/60 
••REDUCE VELOCITY BY '•X < 2) * 3600 
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CHASE (CONTINUED) 



580 PRINT ••D0 YOU WISH T0 TRACK THE SHIPS**! 

590 INPUT AS 

600 IF AS=**YES**THEN 620 

610 ST0P 

620 PRINT ••TRACKING INTERVAL^^l 
630 INPUT J 
640 PRINT 
650 PRINT 

660 PRINT ••TRACKING DATA** 

670 PRINT - 

680 PRINT 

69 0 PRINT **ELAPSED*V**RELATI VE**#**RELATI VE"#**KELATIVE'*#'*TIME T0** 
700 PRI NT *'TI M E*'# "AN GL E**# ''DI STAN CE*** "VEL0 CI TY **# **REN DEZ V0 US*' 
710 PRINT ** **#*• ••# 

720 F0R T = 0 T0 INT(T(2)/60) STEP J 
730 PRINT t# 

740 LET W(3) = T* 60* VC 3) + W( 1 ) 

750 LET W< 4)=T*60*V( 4) + W(2) 

760 LET M = ABS(U(3)-W(/i)) 

770 IF M < 3.14159265 THEN 790 

780 LET M = ABS(M-6. 28318530) 

790 PRINT M* 180/3. 1 41 59265# 

800 LET D = SQR( At 2+Pt 2-2*A*P*C0SCM) ) 

810 PRINT D# 

820 LET V( 5) = ABS( V( 1 ) - V( 2) ) 

830 PRINT V( 5)*3600# 

840 PRINT -T+( T( 1 ) + TC2) )/60 

850 IF T = T(2)/60 THEN 890 

860 NEXT T 

870 LET T = T(2)/60 

88 0 G0T0 730 

89 0 PRINT 

900 PRINT ••EXECUTE INITIAL RETR0FIRE^^ 
910 PRINT 

920 LET L = S*( 1-Et 2) 

930 LET 8 = S*SQR(L/S) 

940 LET W( 1) = W( 3) 

950 LET WC 2) = W( 4) 

960 IF W(l) < 6.28 318 53 THEN 990 

970 LET U(l) = W( l)-6. 28318530 

98 0 G0T0 9 60 

990 IF W(2)< 6.2831853 THEN 1015 
1000 LET W(2) = W(2)-6. 28318530 
1010 G0T0 99 0 
1015LET J=5 

1020 F0R T=( J + T( l)/60) T0 T( 1 ) / 30 STEP J 

1022 LETQ=Q+ 1 

1030 LET K = 0 

1040 LET 0 = 6.2831853 

1050 LET U = 0 
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CHASE (CONTINUED) 



1060 LET Z = <0+U)/2 

1070 LET N<1) = Lt2*SQR<S*62747)/<60*B*62747) 

1080 LET N<2) = E*SIN<Z)/<<Et2-l)*< 1 + E*C0S<Z))) 

1090 LET N<3) = 2/ SQR< 1 - Et 2) t 3 

1100 LET N<4> = ATN< SQR< 1-Et 2)*TAN<Z/2) ) 

1110 LET N s N< 1)*N< 2)+N< 1)*N< 3)*N< 4) 

1113 IF Z<=3. 1 41 59265 THEN1120 

11 15 LET N=N + T< l)/30 

1120 LET K = K^-1 

1130 IF ABSCN-TX.Ol THEN 119 5 

1140 IF K= 1 5 THEN 1 19 5 

11 50 I F N > T THEN 1 18 0 

U60 LET U = Z 

1170 G0T0 1060 

1180 LET 0 = Z 

1190 G0T0 1060 

1195 IF T<T(l)/30 THEN 1200 

1197 LET T< 3) = T< l)/60+T( 2)/60 

1198 G0 T0 120 5 

1200 LET T< 3) = Q*J+T( 2)/60 
1205 PRINT T< 3)^ 

1210 LET W<2) = J*60*V( 4) + W< 2) 

1220 LET W=Z-3. 14159265 +W< 1) 

1230 LET M=ABS< W-W< 2) ) 

1280 PRINT M* 180/3. 1 41 59265^ 

1290 LET R = L/C 1 + E*C0S<Z) ) 

1300 LET D= SQRC Rt 2+Pt 2-2*R*P*C0SCM) ) 

1310 PRINT 

1315 IF T<TCl)/30 THEN 1320 

1316 LET VC 1)= SQRC C 2*S*62747-R*62747)/C S*R) ) -XC 

1317 60 T0 1330 

1320 LET VC 1) = SQRC C 2*S*62747-R*62747)/C S*R) ) 
1330 LET VC 5) = ABSC VC 1) - VC 2) ) 
1340 PRINT VC 5)*3600^ 

1350 LET TC4)==TCl)/60 +TC2)/60 -TC3) 

1351 PRINT TC 4) 

1360 IF T=TCl)/30 THEN 1420 
1370 NEXT T 
1380 PRINT 

1390 PRINT "EXECUTE FINAL RETR0FIRE** 
1400 PRINT 

1405 LET T=TC l)/30 

1406 LET J=TC4) 

1407 LET R=P 
1410 G0 T0 1022 
1420 ST0P 

1430 END 
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ORBITAL TRANSFER 



REF: XFR 



The problem of "where are we?" besets every pilot of a spaceship once 
he has fired his ship's motors. The problem in its most simple form is finding 
the new orbit attained from the old orbit by introducing an energy vector re- 
presented in this program by the coplanar thrust of the motors. 

The program despite its complexity does have simple inputs and outputs. 
The inputs for the old orbit are in response to interrogation. 

1, Eccentricity of the old orbit. 

2, Length of the Semi-Latus Rectum in nautical miles. 

3, Position of the ship in degrees from perigee axis. 

The program then invites the student to fire the motors which he does by answering 
further interrogation . 

1. Angle at which motors are fired (in degrees), 

2. Length of time motors are to be fired (in seconds). 

3. Desired interval between interim calculations (in seconds). 
The output from the program is quite prompt, 

1, Eccentricity of new orbit. 

2, Length of new Semi-Latus Rectum in nautical miles, 

3, Length of new radius vector in nautical miles, 

4, New radius vector angle with respect to the new perigee axis, 

5, Direction and amount (in degrees) of orbital rotation as measured 

between the two perigee axes , 

Obviously, this program can be sequentially staged to move the ship about quite 
freely. Attempting to move under power from one polar plot to another is sub- 
stantially challenging, particularly if velocities are to be matched. Repeated use 
of this program demonstrates that axial shifts are much larger for low eccentricity 
orbits than for ones of higher eccentricity given the same amount of disturbance 
through added thrust. 
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The general factors that must be taken into consideration are presented 
in crude and exaggerated form in the following figures . 
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ERIC 



An apology Is extended for the profusion of factors in Figures 1 through 3, 
but each has its particular place in orbital navigation. Before proceeding to 
a fuller discussion of the program a small development of the use of the 
inclination angle is in order. Certain liberties have been taken in its deri- 
vation which do not appear in the mathematical text and could be puzzling to 
the serious student. 

By assuming that each moment the motors fire is their last moment, and 
consequently at each moment the ship is in a calculable orbit, these inter- 
mediate orbits may be continuously calculated. The program uses an interval 
of 5 seconds but is adjustable to any figure one may desire. Obviously the 
smaller the interval the more accurate the calculations become. In an orbit 
of thousands if not tens of thousands of seconds in period an interval of 5 
seconds closely approximates the actual value. Assuming the 5 second (or 
less) interval is a sufficiently close approximation of the actual value, the 
calculations for instantaneous angle of orbital inclination become rather straight- 



Since we are substituting a small t for dt the angle swept out by the radius 
vector for small t is arbitrarily held to approximate that angle swept out for dt, 
i.e. , d(9 + 9*). The radius vector, r, increases by dr across angle d(9 + 9*) 
and the differential arc length approximated by its tangent has length rd(9 + 9*). 
Assuming d9 then the "new" radius vector crosses the tangent to the curve at 
right angles. By taking the angle vertical to 0 the construction becomes clearer 
as noted on the following page. 



forward. 




Figure 4. 



Angle of Inclination 
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Figure 5. \ 



Tan 0 = dr from use of 

rd(9 + 9*) 
vertical angles 

(NOTE: Error of small t in place of dt is proportional to the difference between 
the arc length of the orbital path and the length of the equivalent tangent seg- 
ment.) 

If r = L then dr = e sin(9 + 9*) by the 

1 + e cos (9 + 9*) rd(9 + 9*) 1 + e cos (9 + 9*) 

use of the simple differential formulas. Since the angle of inclination describes 
the orbit all that remains is the task of fitting all the data and formulas together. 

Since the companion program was to be as realistic as possible, an actual 
spaceship was modeled within the program. The ship has a certain mass, thrust, 
and fuel usage compatible with current chemically propelled rockets. These ship 
parameters can be varied for any number of reasons . 

The ship is at P^: (rQ, 9o) in an originating orbit with known constants of: 

To = ^o (El) 

1 + e^cos 9o 

At this point the ship has a velocity, v^, determined by 

vo^ = Gm (2 - (1 - epZ) ) (E2) 

where Gm is Earth gravitational constant. 




Furthermore, the angle made by the tangent to the orbit and the perpendicular 
to the radius vector at Pq is the angle of inclination and is 



Tanjefo = ^o sin 9o (3) 



1 + eQ cos 9o 



The ship is rotated to the burn angle (measured from the perpendicular to the 
radius vector to the axis through motors) and the motors are fireH. At the end 
of (m) seconds the ship has added perpendicular and radial distance vectors due 
to this thrust. These distance vectors are: 

= Im^ a" cos z (E4) 

and 

= im^Ca sin z) (E5) 
The supporting equations for these vectors are: 

g^ = Gm (E6) 

a • = T^ ' . (E7) 

(w - m . k/2) 

h = mvo sin 0q + jm^ (3" sin z) (E8) 

gh = Gm/(r + h)^ (E9) 



and 



where 



¥ = (go + gh)/2 (ElO) 

t is the ship thrust in pounds; 

w is the ship weight in pounds; 

k is the pounds of fuel burned per second; 

g^ is gravitational acceleration at Pq ; 
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gj^ is gravitational acceleration at height h; 

g is mean gravitational acceleration during time m; 

and 

a" is mean acceleration from ship during time in. 

To keep these equations at a secondary school level of complexity certain 
simplifying assumptions have been made. Evidently g and a are not true 
means but first approximations from their respective infinite series. Also, 
for m, a time interval of m ^ 5 should be used. The results of (E3) are 
based on the tangent to a curve closely approximating the curve for small 
distances; i.e. , the differential arc length is extended. 



At the end of m seconds of firing the ship is in a new orbit which may be 
computed in terms of the old orbital parameters. 

For the point of entry P(r,9) on the new orbit, r. is computed as 

r2 = (xi +X2)2 +(ro + yi +y2)^ (Ell) 

X2 = mvQ cos 0Q (El2) 

y2 = mvQ sin 0q (E13) 

The radius vector has now moved forward angle (w) from its old position. 
Further orbital relationships are: 

9 = w + 9o (E14) 

Tan w = (xj + X2)/r (El5) 

m^v^ = (xi + X2)2 + (yi + vz)^ (E16) 

Tan jar = (yi +y2)/(xi + X2) (E17) 

The original reference line which was the perigee axis of the old orbit 
has rotated to reflect the perigee axis of the new orbit; i.e. , the ship's 
reference line has rotated through an angle of 9*. 

The new orbit through point P is defined by: 



r = 



L (E18) 



1 + e cos (9 + 9*) 
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v2 = Gm (2 - ( 1 - e2) ) (E19) 
r L 

Tan 0 = e sin (9 + 9*) 

1 + e cos (9 + 9*) (E20) 



The simultaneous solution of (E18-20), a tedious exercise in analytic 
geometry, yields: 

e2 - 1 = rv2 cos ^ 0 (rv2 - 2Gm) (E21) 
Gm2 

L = r 2v2 cos2 0 (E22) 
Gm 



sin 9* = -c sin 9 



+ cos 9 ^l 1 - c2 (E23) 



where 

c = (L - r)/er = cos (9 + 9*) (E24) 



A sufficient condition for -1 £ c ^ 1 is that P £ r 6 A and of the two 
solutions to (E23) the sum is used in the first two quadrants while the 
difference is used in the last two quadrants. 
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XFK 



100 


PHINT "THIS PK3GKAM IS DI RECTLV KEYED 13 THE 


T;3PIC 0UTLINE" 


110 


PKINT •••0HB1TAL TRANSFEK' AND THE NUMPEKS IN 


THE 


REM C9LUMN" 


120 


PKINT "^F THE LIST 0F THIS PFx0GKAM REFER Ta THE 


EQUATI3N NUMBERS 


130 


PKINT"IN THE TEXT." 






150 


PKINT 






160 


LET 


T = 75000 'THRUST 






170 


LET 


G = 62747 'GM 






18 0 


LET 


W = 50000 '^'EI GHT 






190 


LET 


K = 100 'FUEL USE 






200 


LET 


I = 180/3.14159265 'RAD/DEG 






210 


LET 


J = 3.14159265/180 'DEG/RAD 






2'^o 


PRINT "ECCEiNTRI CITY"; 






230 


INPUT E 






240 


PRINT "SEMILR"; 'N-M. 






250 


INPUT L 






260 


PKINT "P3SITI0N ANGLE"; 'DEC 






270 


INPUT P 






2«0 


LEI 


P = J * P 






290 


PRINT "BURN ANGLE"; 'DEG 






300 


INPUT Z 






310 


LET 


Z = J * iil 






320 


PRINT "BURN TIME"; 'SEC 






330 


INPUT T( 1 ) 






340 


PKINT "INCrxE.MENT"; 'SEC 






350 


INPUT M 






360 




N = M 13 T( 1 ) STEP M 






370 


LET 


K = L/( 1 + F:*C0S(F)) 




'C El ) 


39 0 


LET 


V = SuR(C) ^^ i:.L)K( ( 2/K) - ( ( 1 - ET 2) /L) ) 




'C E2) 


410 


LE'I 


U = AINC E*blN(P)/( 1 + E*C;3S(P) ) ) 




•CF3) 


430 


LET 


G( 1 ) = G/ ht 2 






450 


LEI 


A = T*G( 1 )/( U-H*K/2) 




•C E7) 


470 


LEI 


H = V* bINC U) + . 5*MT 2*A* SIN( Z ) 


•C EH) 


4y 0 


LEI 


G( 2) = G/ ( R+H) T 2 






510 


LET 


G( 3) = ( G( 1 ) +G( 2) )/ 2 




•C ElO) 


520 


LEI 


A = 1 * G( 3)/( l'.-.M*K/ 2) 




•CE7) 


540 


LET 


X( 1 ) = . 5*MT 2* A* S( Z) 




• C E4) 


5^0 


LE i 


Y( 1 ) = . 5+MT 2*A*SIN(Z) 




•C F5) 


5c50 


LET 


AC 2) = t>^*V*C'*^ S( L) 




•C F 1 2) 


6U0 


LEI 


r ( 2) = V*bIN ( U) 




•C E 1 3) 


620 


LEI 


PC 1) = ATNCCXC 1)+XC2))/K) 




•C E 1 5) 


640 


LEI 


K = SURC CXC 1 ) +XC 2) ) T 2+ C k+YC 1 ) + YC 2) )T ir!) 




•C Fl 1 ) 


660 


LET 


P = p( 1 ) +p 




•C El 4) 


6c50 


LEI 


b = bJRC CX C 1) +X C 2) ) T 2+ C YC n 4^ Y C 2) ) t ?3 




•C F 1 6) 


6^0 


LET 


V = b/.^i 






710 


LET 


U = ATM C C Y C 1 ) +YC 2) ) / CX C 1 ) +XC 2) ) ) 




'C E 1 7) 


730 


LFl 


E = bUhC C h*\/T 2*CJbC U) t 2+ C K* Vt 2- 2+ G) + Gt 2) / Gt 2) 




7b0 


LE i 


L = C C f^★\/★CJ bC U) ) t 2) /G 




•C F22) 


770 


Lf- I 


C = CL-K)/Cb"*h) 




•C F?'4) 


7^^ U 


I F J 


^ < 6. 2H 3 1.-^5 30 THEN c^OO 






'NO 


LET 


f' = p-6. 2n3ldb30 
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XFK (CONTINUED) 

BOO IF P > 3.14159625 THEN 8 40 

aiO LET S= (-C*SIN(P) + C0S(P)*SQR( 1-Ct 2) )/2 

81 1 LET S=2*S 

820 LET S = ATNC S/SQK( 1-St 2) ) 
830 G0T0 9 30 

840 LET S = (-C*SIN(P)-C0S(P)*i>QK( 1-Ct 2) )/2 

841 LET S=2*S 

850 LET S = ATNC S/SQKC 1-St 2) ) 
860 IF Y( 1) < 0 THEN 89 0 
870 LET P = P+S 
875 LET S( 1) = S( 1) - S 

88 0 G0T0 9 70 
B90 LET P = P-S 

89 5 LET S( 1) = S( 1) + S 

900 LET S = ATNC DC 2)/SQK( 1-D( 2) T 2) ) 'AXIAL SHIFT 

9 10 G0T0 9 70 

920 

930 I F Y( 1 ) < 0 THEN 9 60 

940 LET P = P-S 

945 LET S( 1) = S( 1) + S 

950 60 T0 9 70 

960 LET P = P+S 

965 LET S( 1) = S( 1)-S 

970 LET W = W-M + K/2 

980 IF P < 6.28318530 THEN 1000 

1000 LET W = U-M+K/2 

1020 NEXT N 

1080 PRINT 

1090 PRINT "ECCENTRI CI TY"^ "SEMI LR'S "RAD V"* "RAD AN GLE"* "AX I AL SHIFT 
1100 PRINT E^L^ R^P*I ^ S( 1)*I 
1500 END 



~~ - -ni.B^ ^2 
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RUN 

XFR 03/18/70 09: 32 

THIS PROGRAM IS DIRECTLY KEYED T0 THE T0PIC OUTLINE 
•ORBITAL TRANSFER' AND THE NUMBERS IN THE HEM COLUMN 
0F THE LIST OF THIS PROGRAM REFER T0 THE EQUATION NUMBERS 
IN THE TEXT. 

ECCENTRICITY? .01 
SEMILR? 48 0 0 
POSITION ANGLE? 90 
BURN ANGLE? 30 
BURN TIME? 20 
INCREMENT? 5 

ECCENTRICITY SEMI LR RAD V RAD ANGLE AXIAL SHIFT 

7.49 169 E-3 4896.69 4924.29 184. 512 -93.667 

TIME: 0 . 498 SEC. 
READY 
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